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ABSTRACT 

We present an analysis of the dipole straylight contamination (DSC) for 
spinning space-missions designed to measure CMB anisotropics. Although this 
work is mainly devoted to the Planck project, it is relatively general and 
allows to focus on the most relevant DSC implications. We first study a simple 
analytical model for the DSC in which the pointing direction of the main 
spillover can be assumed parallel or not to the spacecraft spin axis direction 
and compute the time ordered data and map. The map is then analysed paying 
particular attention to the DSC of the low multipole coefficients of the map. 
Through dedicated numerical simulations we verify the analytical results and 
extend the analysis to higher multipoles and to more complex (and realistic) 
cases by relaxing some of the simple assumptions adopted in the analytical 
approach. We find that the systematic effect averages out in an even number of 
surveys, except for a contamination of the dipole itself that survives when spin 
axis and spillover directions are not parallel and for a contamination of the 
other multipoles in the case of complex scanning strategies. In particular, the 
observed quadrupole can be affected by the DSC in an odd number of surveys 
or in the presence of survey uncompleteness or over-completeness. Various 
aspects relevant in CMB space projects (such as implications for calibration, 
impact on polarization measurements, accuracy requirement in the far beam 
knowledge for data analysis applications, scanning strategy dependence) are 
discussed. 

Key words: Cosmology: cosmic microwave background - space vehicles: 
instruments - methods: data analysis. 



1 INTRODUCTION 

The observed pattern of the cosmic microwave back- 
ground (CMB) anisotropics is dominated by the dipole 
signal originating from the motion of the Sun in the 
CMB rest frame. This signal is two orders of magni- 
tude larger than the typical cosmological effect. While 
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the dipole subtraction is of primary importance to map 
CMB anisotropics, little is known about the stray- 
light contamination induced by this dipole pattern. The 
straylight contamination consists in unwanted radia- 
tion entering the beam at large angles from the an- 
tenna boresight direction. The most relevant contribu- 
tion is expected from a particular far beam extended 
region, the so-called main spillover, typically located 
at ~ 40° — 90° from the main optical axis. In the 
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case of the Planck ^ telescope, its pea k response is 
at ~ 90° from the telescope optical axis iSandri et alJ 
l2004l ^ and points far away from the Sun along direc- 
tions not far from the spacecraft spin axis. It is al- 
ready known that the straylight contamination of other 
non-cosmological effects, such as Galactic e missions 
iBurigana et al.l2001al . l2004l:ISandri et al.l2004) and in- 



ner Solar System bodies teuriga na et alj[200d. [2001^) . 
plays an im portant role in design i ng space-missions 
hke Planck jMandolesi et alJll998l 120021: iPueet et alJ 



ll998l : lLamarre et alJl2002l:lTaubeill200lLl2006ll . It~ 

trol and removal represents also a relevant aspect of 
WMAP ^ design a nd data analysis iBarnes et aljl20o3 : 
IJarosik et a'i]l2006ll . 

The aim of this paper is to study the dipole stray- 
light contamination (DSC) with analytical and numeri- 
cal methods. We start building a model sufficiently sim- 
ple to be treated analytically that allows to understand 
to first order how the straylight features affect the re- 
coverd CMB angular power spectrum in a way largely 
independent of the specific optical and scanning strat- 
egy details. The analysis is first performed considering 
the direction of pointing of the main spillover centre par- 
allel or not to the spin axis but keeping the the main 
spillover centre direction in the plane defined by the 
telescope and spin axis directions. In this work, we call 
a the angle between the main spillover centre direction 
and the spin axis direction. We then extend the analysis 
allowing for displacements of the main spillover centre 
direction from the plane defined by the telescope axis 
and the spin axis. The analytical results are presented 
perturbatively in a. 

We shall work in the fixed frame with the satellite 
with axes pointing fixed (far away) stars. In this frame 
the vector associated to the dipole is constant while the 
main spillover direction is not (it rotates of 2ti in 1 year). 
We consider the dipole for the motion of the Sun with 
respect to the rest frame of the CMB and we neglect, 
for simplicity, small deviations due to the motion of the 
Earth around the Sun. 

We implement dedicated numerical simulations in 
order to check the analytical results and extend the 
analysis to more complex (and realistic) cases by re- 
laxing some of the simple assumptions adopted in the 
analytical approach. In particular, simulations allow to 
consider non-small values of the angle a and complicate 
scanning strategies introducing the effect associated to 
the survey uncompleteness (or over-completeness) and 
extending the analysis to higher multipoles, I. 

Finally we discuss how the DSC may affect the 
multipole pattern, with particular emphasis paid to 
the possible conne ction with the small s ignal observed 
by COBE/DMR (Hins haw et all Il99d) and WMAP 



http://www.rssd.esa.int/planck 
^ http:/ /lambda. gsfc.nasa.gov/product/map/currcnt/ 



iBennett et alJl2003al: ISoergel et a'i]l2003l l200dl at low 

t 

The article is organized as follows. In Section|21the 
convolution of the dipole and the straylight beam is pre- 
sented. The analytical model for the beam response in 
the main spillover region is presented in Section|^where 
the map is computed with particular care to low multi- 
poles (up to ^ = 4). In Section 2] the power spectrum is 
discussed and the DSC of the quadrupole {i — 2) is anal- 
ysed. In Section 1^ numerical simulations are presented 
and discussed. In Section |S| we extend the analysis to 
arbitrary directions of the main spillover. The main im- 
plications for spinning space missions dedicated to the 
CMB anisotropy are discussed in Section|7| In Section|S| 
a statistical analysis for the amplitude of the quadrupole 
is presented. Finally, our main conclusions are drawn in 
Section H 



2 THE DIPOLE AND THE STRAYLIGHT 
BEAM 

We start considering the convolution / of the dipole 
with the main spillover: 



dnTi,r^Yi"'ie,if)BsL{e,ip), 



(2.1) 



where rffi is the element of solid angle, dfl — d6 sin 6 dip 
with the colatitude 6 £ [0, n] and the longitude ip £ 
[0, 27r [, the sum on m over —1, 0, 1 is understood, Tim 
are the coefficients of the expansion of the dipole ^ 
on the basis of the spherical harmonics Yl^{6,ip), and 
Bsl{0,ip) is the beam response representing the shape 
of the main spillover in the (9, (/p)-plane. In this nota- 
tion BsL is normalized to the whole beam integrated 
response, dominated by the contribution in the main 
beam 



dnS : 



dQ.B ~ 2-Kal 



(2.2) 



main beam 



where at, = FWHM/V8 In 2, with FWHM " represent- 
ing the main beam angular resolution. 

The convolution / can be rewritten in the following 

way: 



Tio / dO dip smO cos 6 BsL{0,ip)- 



\f2 j de dip sm e Re [Tiie^'^] BsL{0,ip) 



(2.3) 



where Re [...] is the real part of [...]. In order to obtain 
equation H2.3^ it has been used that Ti-i = — Tn, where 

^ We use the symbol Tf^. because we want to make clear 
that the dimensionality is given by a temperature (e.g. in 
mK). 

4 FWHM = Full Width Half Maximum. 
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the symbol * means complex conjugation. For sake of 
completeness, we write the spherical harmonics for £ = 1 
l)Sakuraill985l) : 



(2.4) 



e^"^ sin t 



(2.5) 



Although it may not be clear from the notation of equa- 
tion II2.3II . notice that / is a function of the geometric 
features of the shape of the main spillover in the (61, ip)- 
plane. 



3 ANALYTICAL MODELS 

Equation 1)2. 8|l is general and exact. Any specific ap- 
proximation of the beam response Bsl inspired by ex- 
periments will introduce a certain degree of uncertainty 
with respect to the real case. We introduce two different 
parametrizations for Bsl'- we first consider the Top Hat 
approximation, analysed in Subsection l3.1l and then the 
Gaussian approximation, studied in Subsection 13.21 We 
shall see in the following that the two descriptions are 
quite similar and are formally equivalent when the main 
spillover area is sufficiently small (see Subsection 13.31 1. 
For this reason we concentrate on the Top Hat case for 
the analytical estimate of the DSC. 

3.1 Top hat approximation 



Our first approximation iGruppuso. Burigana fc Finellil 
12005 ') for Bsl is the Top Hat: 

BsL{e,ip) = /sL A(e, e^s-'Ao<, 6»™, + Ae>) X 

A( (p, ipms — A,p<, tfi^s + A^>) (3.1) 
with A(a,6, c) = S(a - 6) - S{a-c), 

where fsL is a constant (that is related to the ratio be- 
tween the power entering the spillover and the power 
entering the main beam and it is a number much less 
than 1 - it will be estimated in Subsection 14. H and 
Six) is the step function (or Heaviside function) that 
takes the value 1 for j: > and the value otherwise. 
Equation 13.111 is just an asymmetrical rectangular box, 
in the (6', i/3)-plane, centred around the point {6ms, ^Pms) 
and with sides of length Ae>-|-Aa< and A^>-|-A^<. No- 
tice that the point (9ms, ifims) is the direction of pointing 
of the main spillover centre. This approximation leads 
to: 



dQBsL = fsL 



d6 sin 9 



dip. (3.2) 



Considering that 



d9 sin9 cos 9 = - sin (5 -|- 2Ae<) 



sin {5 + 2 9„ 



d9 sin 



Ae< H cos{S + 2 9„ 



sin {5 + 2Ag)/2, 

dtp cos if — 2 cos (^ms 

< 

sin (^A^< + I) , 
dip sin ip — 2 sin I ipms - 



sin (^A^< + 



with e and 5 implicitely defined by Ae> = Ae< -I- 5 and 
A^> = A^< -f e, we obtain the final expression for the 
convolution: 



sin (5 + 2Ae) sin (5-^2 9ms) (2A^ + e) 



. , " 1 A 5 „ ^ , sin (5 4- 2 Ae) 
4W — Ae + - - cos (<5 + 2 6»™,) ^— ^ 




(A. + I), 



(3,3) 



where the label stands for Top Hat. Here we have 
made the notation lighter setting Ae< = Ae and A^5< = 
A^. If the box is centred (i.e. 5 = Q and e = 0) and the 
direction of the main spillover coincides with the spin 
ax is [i.e. (9ms,iPms) = (7r/2,(/7s)1 then the convolution 
is (|Gruppuso. Burigana fc Finelliil2004^ 



fsL 



Ag + 



sin2A 



lie 



(3.4) 



where we have made explicit the dependence on time 
for ips. Notice that the term proportional to Tio drops 
out in this simple case. 

3.2 Gaussian approximation 

The Gaussian approximation can be described rescaling 
the dipole signal as 



TlmS 



(3.5) 



where the factor e is the Gaussian window function, 
^_f(f+i)^2/2^ specified at ^ = 1 and a is the FWHM for 
the main spillover divided by VS In 2, and adopting as 
main spillover response 

J{9,9ms) 



Bsl{9,ip) = b5{ip,ipms)- 



sm{9) 



(3.6) 



where b (as fsL before) is related to the ratio between 
the power entering the main spillover and the power en- 
tering the main beam. In fact equation (13.611 is a pencil 
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beam while the above rescahng of Tirn mimics the convo- 
lution of the dipole with a Gaussian beam. In this way it 
is easy to compute the following Gaussian convolution 



/G^ be-^ /Ax 

V Stt 



[TwV2cose^s - 2Re [Tiie"^'"^] smdms] ■ (3.7) 



3.3 Link between top hat and Gaussian 
formalism 

We shall see in Subsection 13 .41 that we are interested in 
the following directions for the mainspillover 



cos 9,nb a + O (a 



= ± ^ =F sin 6i„6 a + O (a^) 



(3.8) 
(3.9) 



Replacing equations 13.811 and 13.91 in equation 13.31 
and considering a symmetrical top hat (i.e. e = 5 = 
and — Ae = A), one obtains the following expres- 
sion 



sin 2 A \ 
A H — ) sm A X 



[T,,e-] + Asin(2A) cosga 
^ ^ \/2(A-h 2ii^)sinA 

-Re [Tiie'"^] sin^ o I , (3.10) 



where only the linear term in a is kept. Inserting equa- 
tions 13.81 and 13.91 1 in equation 13.71 and expanding in 
a again to the linear order one finds 



(3.11) 



|±Jm [Tne'*'] + ^^^^lo - Re [Tuc"^] sin^aj 



Equation 13.101 and equation 13.111 are very similar and 
become the same expression for small A, if we define b 
as follows 




be-" =4fsLA^ 



(3.12) 



This is an intuitive result: up to normalizations, a pencil 
beam and a top hat beam are not distinguishable when 
the (angular) area for the top hat is sufficiently small. 
Given this result, in the following analytical treatment 
we will focus, for sake of compactness, only on the sym- 
metrical top hat case adopting 



e - 
At 



S = 0, 

= Ac. : 



(3.13) 
(3.14) 



in the numerical estimates. Instead, we will widely con- 
sider the Gaussian case in the section dedicated to the 
numerical simulations. 




Figure 1. Cone drawn by the direction of the main spillover 
(O-MS) during the rotation of the main beam (O-MB). 0^i, 
is the angle between z axis and O-MB, tps is the angle in x-y 
plane between x axis and spin axis (O-S) and a is the angle 
between O-S and O-MS. 



3.4 Relation between main beam and main 
spillover 

During the rotation of the main beam, the main 
spillover, if not lying on the spin axis, draws a cone, 
as depicted in Figure This means that the direction 
of the main spillover and of the main beam are related 
by 



dms ~ — — arctan(tana cos^mt) 
fms = </5s ± arctan (tan a sin 6m.b 



(3.15) 
(3.16) 



where Omt is the colatitude of the main beam, a is the 
angle at the vertex of the cone and ips is the longitude of 
the spin axis (see also caption of Figure . The upper 
(lower) sign of equation 13.161 has to be taken into ac- 
count when the main beam rotates from North (South) 
to South (North). 

In order to make realistic this simple model, we 
have to check that during the rotation, the solid angle, 
subtended by the main spillover (Qms), is constant. A 
simple calculation gives 

= 4A sin A sin 61^^ = 4A smA ^ ^^^^^ 

Vl + tan a cos^ 9mb 

Of course equation l|3.17|l is not constant because 6mb 
is a function which depends on time. However, we have 



: 4A sin A 



-cos 



0„6 +0 (a*) 



(3.18) 



and therefore ^Ims is constant in time to linear order 
in a. This means that the analytical calculation will be 
done to linear order in a. Then equations 13.151 and 
13.161 will be Taylor expanded for small a, obtaining 



"ms — 2 



± sin ( 



(3.19) 
(3.20) 
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Figure 2. Sketch of the nominal scanning during the first 
survey. 

3.5 Building the map 

The total signal received by the satellite receiver, is the 
sum of the two following terms: 



T{e,ip)^TMB{O,ip)+IsLi0,'p). 



(3.21) 



here Tmb is the signal entering the main beam (where 
the dipole has been subtracted away) whereas Isl is the 
signal due to the dipole entering the main spillover. We 
take into account for simplicity a simple scanning strat- 
egy as that considered in the Planck project as nom - 
inal scanning strategy (NSS) llDupac fc Taubeil [20051') . 
i.e. with the spacecraft spin axis always in the antiso- 
lar direction or, in practice, always on the ecliptic plane 
as valid when the small effect of the spacecraft orbit 
around the Lagrangian point L2 of the Sun-Earth sys- 
tem is neglected. 

In this case, for the first survey Isl is given by 



I' 



- a cos 6,(p + ^ — a sin i 
for < If < n 



[t: — ct cos 



?, (/p — § + sin 9^ 



(3.22) 



for TT < if < 2tt , 



whereas, for the second survey it is 



— a cos tl, ip - 



+ a sin ^ 



for < < TV 



I^" (f - a cos 6^, <^ + f - a sin f 
for n < Lp < 2n . 



(3.23) 



The shift in the definition of the arguments in Isl (dur- 
ing either the first or the second survey) comes from 
the fact that when the main beam rotates from North 
to South the main spillover is shifted of — 7r/2 plus a 
small correction proportional to a (due to the non- 
perfect alignment of the main spillover with the spin 
axis), while when the main beam rotates from South 
to North the main spillover is shifted of +'k/2 minus 
a small correction proportional to a (still due to the 
non-perfect alignment of the main spillover with the 



spin axis). Moreover note that it has been used that 
~ V'mb— 7r/2 when the main beam rotates from North 
to South while (ps = fmb + 7r/2 when the main beam 
rotates from South to North. See Figure |5| for a sketch 
of the scanning. 

Notice that now {9, ip) is the pointing of the main 
beam (we omitted the label mh to make the notation 
lighter) . 

3.6 Definition of T^^ 

As usual, we expand the signal in spherical harmonics 
T[9, = ^ Te^Ye^{9, <f) , (3.24) 

Im 

which implies 



T,^ = / dflT{9,ifi)Y;^(9,ip). 



(3.25) 



We consider the multipole expansion of the straylight 
contribution 



TfJ^ = / dn Isl{9, >f) y;u9, <p) . 



(3.26) 



To make lighter the notation, we rewrite the convolution 
in a compact form as 

l'^" {9ms;'Pms) = Cl Sm{29ms) — 

(C2 — C3 COs{29ms)) (dl COS (fims — ^2 Sinifms) (3.27) 



with 



Cl = \/3/47r/sLA sin(2A)Tio 
C2 =4y378^/siA 

C3=4y^3/8^fsL sin(2A)/2 
di = sin A Re [Tn] 
d2 = sin A Im [Tn] . 



(3.28) 
(3.29) 
(3.30) 
(3.31) 
(3.32) 



3.7 General results after an odd number of 
surveys 

After one or an odd number of surveys it is possible to 
show that all the TfJ;^ with odd m vanish when a = 0, 
as follows from: 



rpO _ 

J-im — 



d9 sin 9 



TV TT J. 



(3.33) 



where the label stands for q = and where the label 
has been dropped for simplicity. Changing variable 
of integration [ip = ip' — n in the first integration and 
ip = (p' +TT in the second integration of equation I3.33ll 1 
and considering that 



YL(9,p±7t) = {-irYU9,p)., 



(3.34) 
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then equation 13.331 reads 

[i-(-i)'"]r«„ = o. 



(3.35) 



Equation H3.35|l is identically satisfied when m is even 
and demands = when m is odd. 

Moreover, still considering equation H3.33^ . the in- 
tegration over X = cos6' vanishes when Yfm(arccos a;, i^j) 
is odd with respect to x. This happens when i is odd 
and m is even ^, i.e. : 



if ^ is odd and m is even . 



(3.36) 

Equations 13.351 and 13.3611 show that only Tj,'^ with £ 
even and m even, can be non- vanishing. In particular 
Tirn = if is odd, a property that will be widely used 
throughout the article. 

Considering a small but not zero then is possible 
to write Tf^ as follows 



TfiS L rriO 



ot r/„ , (3.37) 

where the label ^ stands for first order in a. It turns out 
that 

1 <.2ir 



dx / di^ 
1 Jo 



2cix 



(c2 + C3 ) (di cos if — d2 simp) \/\ 

Yl^{&tCCQ%X,if>) , 



(3.38) 



where x — cos 9. For symmetry reasons it is possible to 
show that only 



Ten 



1 f2Tr 

dx / d(p2cixY(Q{a.Tccoax,ifi) , (3.39) 
1 Jo 



and 



T/±i ^ J dx J 



1 p2TT 

dx / d(y3(c2 + C3)(di COS — d2 sin (p) 
1 Jo 



•\/ 1 — x'^Yi-^i{a,Tccos x,ip) . 



(3.40) 



can be different from zero when £ is odd. It turns out 
that both equations 13.391 and 13.4U1 vanish except the 
case £ = 1. We calculate 



9 



Tlo - 
and 



Tiii = \l + ca) (±di + idi 



(3.41) 



(3.42) 



This means that the unique term that turns on (linearly 
in a) when a is taken different from zero (and small) 
is the dipole term. All the other terms deviates from 
Of = results at least for . 

^ Y(^(a,Tccos X, ifi) is odd with respect to x also when £ is 
even and m is odd. 



3.8 Surviving systematic effect after an even 
number of surveys 

After two (or an even number of) surveys the average 

f,^ = {T^2+T^"^)/2, (3.43) 

where the labels and "-^^^ stand for first and second 
survey respectively, is different from zero only for £ — 1 
(to the linear order in a) . This term vanishes when a — 
0. This means that the whole effect averages to when 
a = 0. 

In order to show this statement, consider that with 
replacements similar to the ones performed in Subsec- 
tion |0 one obtains 

= alem , (3.44) 

where Iim turns out to be the same integral given in 
equation 13.381 and studied in the Subsection 13.71 No- 
tice that the average over two surveys is equivalent to 
one survey result. 

This demonstrates the statement made at the be- 
ginning of this Subsection ^. 

3.9 The computation of 

In Subsection l3.7l it has been shown that the calculation 
of Tf^ to linear order in a is equivalent to the one 
with a = (i.e. T^^) ff we exclude the case £ = 1. 
Moreover in order to have the possibility of obtaining 
a non- vanishing T^^, we must constrain £ and m to be 
even. 

3.9.1 Tcfo^ 

Specifying £ = and m = we obtain 

Tcfo^ = ^di(c2 + cs)+0 (a^) . (3.45) 

The monopole deviates from the order a = as a^. 

3.9.2 TfJ: 

Specifying £ — 1 and m = we obtain 



J 10 ^^y"3"^i"- 

For £ — 1 and m — ±1 we have 



(3.46) 



(3.47) 



° The fact that DSC has null average after 2 surveys when 
a = can be shown also integrating over the time interval 
At corresponding to 2 surveys. 
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Equations 13.461 and 13.4711 show that only in the sim- 
plified case in which the spin axis is parallel to the di- 
rection of the main spillover (i.e. a = 0), we have that 
the DSC map does not show dipole contribution ^. This 
could have some important consequences on calibrations 
based on the kinematic dipole. It is interesting to note 
that this contribution survives even after two surveys. 
This is the unique contribution that does not average 
to zero (when Q 7^ 0) after an even number of surveys. 

3.9.3 Ti^ 

Specifying 1 = 2 and m = 0, it is easy to obtain 

r#o^ = O {cP) , (3.48) 

implying that it is vanishing at linear order in a. Notice 
that even if symmetry arguments would make possible 
a contribution (i.e. £ and m are both even) at the 0*'' 
order in a, the result is vanishing. As it is already known 
from Subsection 13 . 71 for m = ±1: 

TiA = , 



(3.49) 



where no expansion in a has been performed in order 
to obtain this result. For m = ±2: 



To 




(di ± 2id2) (C2 + C3) 



(3.50) 



where a does not appear to linear order [there are cor- 
rections of order O (q^) ]. 

These results have been computed using the defi- 
nition of spherical harmonics for i = 2, tha t we report 
here for sake of completeness (ISakurailll985^ : 



(3.51) 
(3.52) 




r-(.,.)^y^e--sin^.. 



(3.53) 



Equation 13.501 is the (non- vanishing) contribution 
to the quadrupole due to the dipole entering the stray- 
light. This is one of the main result of this work. 



3.9.4 Ti^ 



Setting = 4, m = and m = ±1, we obtain: 

0. (3.54) 



T^S L rriSL 
40 — J-A+^ 



Notice that even if symmetry arguments do not prevent 
a contribution, for 1 = 4: and m = Q, nevertheless the 



^ We will see in Subsection 15.2. fl that equations 113.461 and 
can be generalized for non-small a substituting a with 

sin a. 



computation gives 0. For m = ±2: 

ifii = -J^{d^± 2id2) 4 (C2 + C3) + O (a") (3.55) 
V ivr 15 ' 



For m = ±3: 
Tfii =0. 
For m = ±4: 

.SL 12 /35 



(3.56) 



n±i = ^ ('^i ± 4id2) (C2 + Ci) + [a') (3.57) 

These results have been computed using the defini- 
tion for of spherical harmonics for ^ = 4, tha t we report 
here for sake of completeness JSakurailll985ll : 



160F 



(3 - 30 cos^ 61 + 35 cos* 6) , (3.58) 



F±^(e,v.) = T^yJe±»^x 

cos 6 sin 6 (—3 + 7 cos^ 



y±^(e,^) = |^|-e±^-x 

sin^ e (-1 + 7cos^ ( 



(3.59) 



(3.60) 



F4±='(6',(^) = =f|J— e±^'^cos0sin^e, (3.61) 

8 V 



(3.62) 



4 IMPACT ON THE ANGULAR POWER 
SPECTRUM 

The quadrupole presents a contamination that does not 
depend on a (to the order we are performing the com- 
putation). In this Section we perform a detailed analysis 
of this effect. 



4.1 Comparison among Cf" with low 

We consider 

JdQBsL _ p 



fsL 



(4.1) 



4 A sin A 4 A sin A ' 
where p is the relative power entering the main spillover 
with respect to the total one (i.e. essentially entering the 
main beam). By the definition of Cf^, 



SL 



^ 2^ + 1 ^ 

m— — 

we compute 



J-Pm. ] J- Pi 



(4.2) 



(4.3) 
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2 2 



cos''AT?o + lf{A)\T,,\^ 



04 — 



where 



( f{A) ( 12,^53 : 



(4.4) 
(4.5) 
. (4.6) 



./(A) = (1 



sin A >2 
— ^ — cos A 



As an example, we choose A = vr/lO and p = 1/100. 
Moreover it is possible to show that 



Im [Til] = sin(^d sin6'd 



Tio =cos0d J^T, 



7?e [Ti: 



- cos Lp^ sin 



(4.7) 
(4.8) 
(4.9) 



where {9d,^Pd) is the direction and T is the amplitude 
of the dipole ® . 

In ecliptic coor dinates, according to WMAP 
llBennett et al.ll2003a^ . {ed,<Pd) = (1.7651 , 2.9941)rad 
and T = 3.346 mK. We then obtain: 



/m[Tii] = 0.69823 mK, 
Tio = -1.32225 mK, 
Re [Til] = 4.69963 mK . 



(4.10) 
(4.11) 
(4.12) 



Now we give, as an axample, some numerical re- 
sults. We have 



C^^ = 5029.9 ^lK^ 
Ci^ = 202.74 AtK^ 
Cf^ = 18.222 ^K^ , 

and for a = 7r/36 

Cf ^ = 11.135 ^K^ , 

while for a = 7r/18 

Cf^ = 44.539 ^K^ . 



(4.13) 
(4.14) 
(4.15) 

(4.16) 

(4.17) 



4.2 The angular power spectrum 

From equations 13.2111 and 13.251 it is clear that: 

rT~i rriS , rriS L 



(4.18) 



* These relations are obtained solving the following set of 
equations T = TimY{^(9d,'Pd), = Ti^Y{^i8d + 7r/2, ipd) 
and = TimYj™(7r/2, ipd + '7t/2). Of course, the solution can 
be verified replacing in T = T{ed,'Pd) = TimY{^(ed, <Pd)- 



where is already defined in equation 13.2611 and 

rpSKY ■ 



Tim = j dflTMB{0, f) Yt*m{0, If) . 

The CMB power spectrum is given by 



2£ 



— y 

+1 ^ 



TlmTlm , 



(4.19) 



(4.20) 



and replacing equation 14.1811 one obtains 

I 

ri 1 \ " \rrSKY , rrSLl* [rrSKY , rrSLl 

— 2£ _|_ / , [J-lm +J-lm\ [J-lm + J- Im. \ 

e 

_ 1 \ ^ I /rrSKY\* rrSKY , 

^ 2ZTT ^ ' " 

f rr-\S L \ rj-iS L , / rj-iS L \ rf~iS Y . / rj-iS 

[-•-tm ) J-tm + \J-tm ) J-lm + \J- l-n 



SL-\ 



_ (^Sify ^ qSL ^ ^SKY-SL 

In particular for the quadrupole we have 
Cl'' = \fI^ {[Re [Tn]]' + 4 [Im [Tn]]^) 



) Tt 
(4.21 

(4.22) 
(4.23) 

where it has been used Re [T2±2] = Re [T2^2\ and 
Im [r2±2] = —Im [T2:f2], with Fsl defined as 



and 



C: 



SKY-SL 



= - -Fsl {Re [Tn] Re [Tj^^^] + 
2/m [rii]Jm [Tai^^]) 



/s L — [A + cos A sin A] sin A 
37r 



V5 r 

' 37r 



1 + cos A- 



sin A 



A 



(4.24) 



From equation 14.2411 . choosing A = tt/IO and consider- 
ing p = 1/100 we have Fsl{p = 1/100) ~ 4.59 x 10"^. 

4.3 DSC of the quadrupole 

We are ready to estimate the order of magnitude of ^ 
and C2 . Since T2V„ are stochastic numbers with 

vanishing mean and standard deviation equal to Ci''^^ , 
we adopt, as an example, the following relation 

SKY-\ 2 ^ , r^Sfrvn 2 



2 Re [Tl"''Y^2Im [rg'"'']- 



(4.25) 



for numerical estimate, where the factor 2 is due to the 
assumption that the real and imaginary part give the 
same contribution ^. We obtain: 



= 202.7 ^K^ 



(4.26) 



^ This choice on the real and imaginary part of T^^^ is 
limited to this Subsection. In Section |S] a statistical analysis 
will be performed. 
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±(386.1 ± 114.7) ^K^ 



(4.27) 



where it has been chosen Cf^^ ~ lO^^K^. The ± in 
equation H4.27|l is due to our ignorance about the rela- 
tive sign of Re [tSz'^^] and Im [TjI,^'^] . 

In Table □ there are Ci^ and C^"^^''^^ for p = 
1/500,1/100,5/100 and Cf^'^ = 500, 1000, 1500 /xK^ 
All the numbers for C^^^~^^ contribution have to be 
understood with a ± in front of them (since we do not 
know the total sign of this contribution). 

We stress two main observations. First, we notice 
that the contribution is smaller than the C|^^~'^^ 
contribution if p is sufficiently small. This is clear be- 
cause ^ is quadratic in p while C|^^^~'^^ is linear. 
Therefore, for p = 1 /500 and p — 1/ 100 we have that the 
leading term is given by C^^^^^^ while for p = 5/100 



we find that ^^"^^ 



is subleading; second, since Cl 



contribution is always positive, because of the previous 
observation, it is clear that Cf can be lowered only 
for the first and second column of Table when the 
relation 14.2511 is chosen. 



4.4 Decrease of the quadrupole 

From equation 14.2111 we can rewrite the observed 
quadrupole C2 as a function of Fsl 



5 5 



c2{Fsl) = cr 

where 

B = Re [Til] Re [ril^^] + 2/m [Tn] Im [T, 



SKYl 
22 



A=[Re [Tii]Y+4 [/m[Tii]] 



(4.28) 

(4.29) 
. (4.30) 

Equation H4.28^ is just a parabolic behaviour in Fsl 
with concavity in the upward direction since the coef- 
ficient of Fgj^ is always positive. The sign of B is not 
known a priori: from equation 14.231 or equation 14.281 
is clear that Cf can decrease only if B > 0. We focus 
now on this case. 

The minimum value C2{Fsl\v) is reached by the 
vertex Fsl\v, Fsl\v = B/A, where 

2 B'^ 



C2{Fsl\v) 



5 A 



(4.31) 



2 CRe[rn] Re [T^/^] + 27m [Tn] Im [T^/^] 
5 



^^^'2 ^ <o,5ci"'^/2], we can 



[i?e[rii]]"+4 [Im[T-n]f 

Considering that \T22"" T G 
rewrite equation 14.311 as 

-.SKY _ y 

'^SKY\2 l/KnSKY 



COS^ 7I 



(4.32) 



C2{FsL\v)^Ci' 

lr2l,^'^|V(5C|''«V2) (then / e [0,1]) and 



where / 

COS 7 = [v ■ u)/v with V being the norm of the vector v, 
v^{Re[Tii],2Im[Tii]) , (4.33) 
and u — (cos a, sin a) with a defined as 



a = arctan (Jm [T#2^^] /Re [ri?,^^] ) 



(4.34) 



0.8- 



0.6- 



0.4- 



0.2- 
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-0 . 5 
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Figure 3. Contour plot for the ratio C2(Fsl\v)/C2 
sus the couple of parameters (7, /). Since we are dealing with 
B > than 7 G [— 7r/2, 7r/2]. Each contour represents the 
value n/10 with n G [1,2, ...,9]. The red contour represent 
1/5, the green contour 1/2 and the blu contour 9/10. 



and therefore ranging in the set [0, 27r]. In Figure |21 we 
show a contour plot for the ratio C2{Fsl\v)/C2^^ ver- 
sus the couple of parameters (7,/). Darker is the grey 
and smaller is the ratio (see also the caption). 

In the same way we rewrite equation l)4.28|l as 



y{Fs 



C2{Fsl) 

(JSKY 



f 



1/2 V COS7 



V5Q 



SKY 
2 



■.FSL + - 



5 Ci'^^ 



Fsl ■ 



(4.35) 



(4.36) 



This means that the first branch of DSC is given by 
7 G [— 7r/2,7r/2] (where also a decreasing of the ob- 
served quadrupole is possible, i.e. B > 0) and the sec- 
ond branch is given by 7 £ [7r/2, 37r/2] (where only an 
increasing of the observed quadrupole is possible, i.e. 
B < 0). In FigureBlwe plot equation for Ci'^^ = 

200, 500, 1000, 1500 fiK^ (red, green, blu and yellow lines 
respectively) with (7,/) ~ (tt/IO, 1/2) (lower curves) 
and with (7,/) = (tt/IO + 7r/2, 1/2) (upper curves). 
Bottom panel of Figure 2] is a zoom of the top panel 
but with p as independent variable instead of Fsl (see 
also the caption). 



5 NUMERICAL SIMULATIONS 

The results presented in the previous sections are based 
on an analjd^ical treatment of the DSC effect. Numerical 
simulations allow to independently check the analytical 
results of the previous Section and to extend the anal- 
ysis to more complex (and realistic) cases by relaxing 
some of the simple assumptions adopted in the previ- 



© 0000 RAS, MNRAS 000, 000-000 



10 C.Burigana, A.Gruppuso and F.Finelli 



C2{SKY) 


P = 


= 1/500 


p = 1/100 


p = 5/100 


C2 


500 




8.1 


202.7 


5068.5 


SL 


500 


54.f 


3 ± 16.2 


273.0 ±81.1 


1364.9 ± 405.6 


SKY-SL 


1000 




8.1 


202.7 


5068.5 


SL 


1000 


77.; 


I ±22.9 


386.1 ± 114.7 


1930.3 ± 573.6 


SKY-SL 


1500 




8.1 


202.7 


5068.5 


SL 


1500 


94.6 ±28.1 


472.8 ± 140.5 


2364.1 ± 702.5 


SKY-SL 



Table 1. All the temperatures are given in pK? . See also the text. 



1 . 4 
1 . 2 
1 

. 8 
. 6 




0.002 0.004 0.006 0.008 0.010.012 



1 . 4 

1 . 2 




0.0020.0040.0060.008 0.01 



Figure 4. Top panel: y = y{FsL) for Cf^'*' = 200 /.tK^ (red 
line), C|^^ = 500AtK2 (green line), C|^^ = 1000 ^K^ 
(blu line) and C|^^ = 1500 pK^ (yellow line). Lower curves 
correspond to (7,/) = (tt/IO, 1/2) while upper ones to 
(7, /) = (vr/lO + 7r/2, 1/2). Bottom panel is a zoom of top 
panel but plotted in terms of p instead of Fsl (see equa- 
tion iO^ I: y = y{p) for C§'^^ = 200 ^tK^ (red line), 
C^KY ^ 500AtK2 (green line), C|^^ = 1000 /xK^ (blu 
line) and C|^^ = 1500 /xK^ (yellow line), where A = tt/IO. 
Lower curves correspond to (7, /) = (tt/IO, 1/2) while upper 
ones to (7, /) = (tt/IO + tt/2, 1/2). See also the text. 

ous sections. In particular, we will investigate through 
dedicated computations some aspects that are difficult 
or impossible to address in analytical way: i) the case of 
a non-small angle a (i.e. of a deviation of some tens of 
degrees from the orthogonality between the directions 
of the main beam and the main spillover); ii) more com- 
plicate scanning strategies including, for instance, pre- 
cessions of the spacecraft spin axis around the nominal 



case of a spin axis direction kept along the antisolar di- 
rection, i.e. in practice on the ecliptic plane. Moreover, 
simulations allow to take into account iii) the effect in- 
troduced by the uncompleteness of one of the two sur- 
veys and iv) to extend the analysis to higher multipoles 
in the cases i) and ii). 

5.1 The numerical code 

We used here an update dedicated version of the code 
implemented and successfully tested in many simula- 
tion works devoted to the study (and reduction) of var- 
ious classes of systematic effects in the context of the 
Planck mission under different scanning strategy as- 
sumptions. It is described in d etail in [BuriEana ct aj 
11199811 and iMaino et al.l il999l) and, in particular re- 
gardi ng the straylight effect, in lBurigana et all ll2001al . 
I2OO4I) ■ where further informations on the relevant ref- 
erence systems adopted in the code can be found. We 
do not consider here the impact of the spacecraft orbit 
(around the Lagrangian point L2 of the Sun-Earth sys- 
tem, as for example in the case of WMAP and Planck) 
because its effects are not relevant in this context. 

In this work, we compute the convolutions between 
the main spillo ver response and the sky dipo le signal 
as described in iBurigana et alJ ll2001al . |20o3): but by 
pixellizing the sky at ~ 1° resolution considering 
spin-axis shifts of ~ 1° every day and 396 samplings 
per scan circle, and by adopting the analytical Gaussian 
description of the main spillover response introduced in 

^'^ The HEAL P lx sc heme I http://healpix.jpr.nasa.gov/ 1 by 
iGorski et alj l2005l) has been adopted in simulations. A 
dipole input map at n^i^^ = 64 has been considered. 
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Figure 5. Maps of DSC in the case of the nominal scanning strategy (NSS) and for a = 0°. The parameters of each case 
arc indicated above each map together with the reference to the panel in Figure 13 where the corresponding APS is displayed. 
Galactic coordinates and Mollweide projection are used. 



Section [3.21 that has been exphcitely implemented in 
this version of the code. 

We consider here as reference case a channel at 
70 GHz, a frequency where the foregrou nd contamina- 
tion at large angular sca les is minimum ijBennett et alJ 
l2003bl: IPage et alJl200el) . and report simulation results 
in terms of antenna temperature, as typical in simu- 
lation activity being it an additive quantity with re- 
spect to the sum of various contributions. For numeri- 
cal estimates, we report here the results referring to the 
case p = 0.01 (i.e. a relative power of 1% entering the 
main spillover) and assume a main spillover FWIfM, 
FWIfMms, of 20°, an angular size amplitude compara- 

This is certainly a simple approximation of the complexity 
of the realistic shapes predicted for the main spillover (and 
in general for the beam response far from the mai n beam) 
by optical simulation codes (see e.g. ISandri et al.l j2004)). 
On the other hand, the details of the main spillover response 
shape depend on the considered optical system. In the con- 
text of the Planck activities, they will be included in a fu- 
ture work. The main contribution of this paper is in fact 
the understanding of the most relevant effects introduced by 
the DSC, common to relatively different optical systems, in 
terms of simple parametrizations. 



ble to those of C MB space experiments dSandri et al.l 
I2OO4I : iBarnes et a l. 2003). 

The adopted resolution allows us to investigate the 
DSC effect on angular scales larger than few degrees, i.e. 
up to multipoles ^ ~ 50, where the most relevant effects 
are expected because of the large scale features con- 
sidered here both for the signal and the main spillover 
response. The main output of the simulation code con- 
sists of time ordered data (TOD) containing the signal 
entering the main spillover for each simulated pointing 
direction. We generate TOD for two complete sky sur- 
veys. The TOD are then coadded to produce all sky 
maps of DSC at a resolution of ~ 2° (risidc = 32 in the 
HEALPix scheme). These maps are analysed in terms of 
angular power spectrum (APS) by using the anafast 
facility of the HEALPix libraries. For simplicity, we con- 
sidered always an angle of 90° between the directions 
of the spacecraft spin axis and of the main beam cen- 
tre (assumed to be aligned with the telescope line of 
sight). This assumption allows to reach the all sky cov- 
erage for the whole set of considered scanning strate- 
gies. Differently, in the absence of spin axis displace- 
ments from the ecliptic plane a small unobserved cir- 
cular patch appears around each of the ecliptic poles 
(fPupac fc Tauber..2005i) . We prefer to avoid the case of 
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Figure 6. Maps of DSC in the case of the nominal scanning strategy (NSS) and for o 7^ 0° and in the case of a cycloidal 
scanning strategy (C) with a semi-ampHtude of 10° and a = 0° or 10° . The parameters of each case are indicated above each 
map together with the reference to the panel in Figure |S] where the corresponding APS is displayed. Galactic coordinates and 
MoUweide projection are used. 



non-complete sky coverage because it may generate a 
certain complication in the APS analysis, introducing a 
mix between the effects from DSC and partial sky cov- 
erage. Of course, this choice allows also to check the 
analytical results in the same work condition. 

Our general results are obviously not significantly 
dependent on the specific choice for FWHMms . In addi- 
tion, they can be easily rescaled to any value of p (signal 
oc p, APS oc p^). 

Finally, we note that the discretization of ~ 1° 
adopted in the present simulations implies a typical nu- 
merical error amplitude significantly less (because of av- 
eraging) than ~ 3 mK X (0.5°/180°) X p ~ 0.1 /iK, 
a value at least one or two order of magnitude smaller 
than those relevant in this analysis, and then fully neg- 
ligible for the present purposes. 

5.2 Numerical results 

We carried out the simulation described in the previous 
section for various scanning strategies and configura- 
tions. 

First, we considered the scanning strategy (the 
NSS) previously adopted in Section|21with the spin axis 



always in the antisolar direction for three choices of the 
angle a: 0°,10° and 45°. The corresponding maps for 
some representative cases are shown in Figures |S] and 
|S] We considered the maps obtained coadding the TOD 
from each of the two surveys separately and from the 
two surveys together. In the case with a = 0° we con- 
sidered also the maps obtained from one whole survey 
plus a second uncomplete (2/3) survey. 

We then considered a cycloidal scanning strategy 
(C) with periods of 6 months (i.e. a number, Nc, of 
complete cycloids per 360 days equal to 2), 2 months 
(i.e. Nc = 6), and 1 month (i.e. Nc — 12). The corre- 
sponding maps for some representative cases are shown 
in Figures |S| and |7| In the cases with Nc — 2, we con- 
sidered a = 0° and 10°; in the cases with Nc = 12, we 
considered q = 0° and 45° ; in the cases with Nc — 6, we 
considered a = 0°. The semi-amplitude of the preces- 
sion cone has been set to 10°. For comparison, we con- 
sidered also a semi-amplitude of 5° in the case a = 0° 
and Nc = 12. Again, we analysed the cases of each of the 
two surveys separately and of the two surveys together. 

We observe that it is preferable to first produce the 
map corresponding to each survey separately and then 
to construct the combined map from the two surveys as 
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Figure 7. Maps of DSC in tlie case of a cycloidal scanning strategy (C) with a semi-amplitude of 10° , a = 0° or 45° and for 
Nc = 12 or 6. The parameters of each case are indicated above each map together with the reference to the panel in Figure Isl 
where the corresponding APS is displayed. Galactic coordinates and MoUweide projection are used. 



a simple average of the two separate maps. In fact, even 
a relatively small difference in the number of hits per 
pixel in the two surveys (generated by a combination of 
pixel shape and effective observational strategy) can in- 
troduce a remarkable spurious unbalance in the average 
when it is done by jointly treating the TOD from the 
two surveys 

The APS obtained analysing the various maps are 
reported in Figure |H| for some representative cases. In 
each of the maps in Figures |3|S||7| we report the corre- 
sponding panel in Figure |H1 

5.2.1 Simple scanning strategies 

The simple view to the maps obtained in the case of the 
NSS confirms the results derived through the analytical 
approach. 

Note also that considering a single survey, the numerical 
error in the TOD quoted above is further dropped by av- 
eraging. This is not strictly true by combining two surveys 
because pixel shape and effective scanning strategy imply a 
not perfect geometrical symmetry between the TOD from 
the two surveys and the presence of (very small) residual 
deviations from the ideal case. 



Figure |K| clearly shows the quadrupole produced by 
the dipole straylight contamination when only a single 
survey is considered. Note that the pattern is symmet- 
rical in the two surveys, as expected. 

Considering the average of two surveys, the result- 
ing map is essentially zero, except for a small ring, cor- 
responding to the edges of the surveys, where it is in 
practice impossible to exactly balance the signal from 
the two surveys as a result of the combination of pixel 
shape and efTective observational strategy. 

If one of the two surveys is not complete (see Fig- 
ure|KJ , a remarkable quadrupole again appears. Figure|H| 
clearly shows this effect in terms of APS. 

Note that, as expected from the analytical treat- 
ment for Q = 0, our computation gives an APS different 
from zero only for the even multipoles. 

Note also the power decreasing for increasing multi- 
poles. At low multipoles, we have checked that the APS 
found numerically agrees with the analytical prescrip- 
tion. 

Differently from the ideal analytical case, the 
combination of two complete surveys implies a non- 
vanishing APS because of small deviations of the ex- 
act balance found analytically. In any case, these effect 
(that is found to decrease by increasing the resolution 



© 0000 RAS, MNRAS 000, 000-000 



14 C.Burigana, A.Gruppuso and F.Finelli 



A 

- A. 


-*A. 

*AA A aa 
□ ° ° ° ° ° = □ A A4'"S>*ft 






qI 


b1 


1 1 1 1 

A 

A *^^***^AAAAAAAAAAAAA 


1 1 1 1 

*A 
A A 

*Aa 

° ° ° ° ° "^^^^^^^-i-AfAAA^.^.^ 


n □ ° ° ° ° ° ° . 
=°° A^ A******* 

, A • * " n / 

^ ^ A U A _ 


AAA + + „ □ ° ° ° A A AAA" 
II* 1 1 1 1 


1 1 1 1 

-% 

*Aa A A . 

^'^^^AAAAAAA^AA^^AA^ 


-\ 

A .. 
^A 


***** 

4 """^ □ □ □ ° 


+A++++trxUflsi&fli 


a J 


Jfl°A 

b3 



10 20 30 40 

Multipole (i) 



10 20 30 40 

Multipole {£) 



10"* 



10"* 



c2 



c3 



10 20 30 40 

Multipole (l) 



Figure 8. APS of the DSC for various configurations. Panels al, a2, a3 refer to the NSS, respectively for a = 0°, 10° and 
45°. The others panels refer to cycloidal scanning strategies. Panels bl and b2 refer to a semi-amplitude of 10° and Nc = 2, 
respectively for a = 0° and 10°. Panel b3 refers to a semi-amplitude of 10° and Nc = 12, and a = 45°. Panels cl, c2, c3 refer to 
a = 0° for a semi-amplitude of 10° and Nc = 12 (panel cl), a semi-amplitude of 10° and Nc = 6 (panel c3), a semi-amplitude 
of 5° and Nc = 12 (panel c2). In all panels triangles and crosses refer to a single survey and squares to the combination of 
the two surveys. In panel al, asterisks refer to the case of one complete survey combined with an uncomplete (2/3) survey. See 
also Figures and|7| Note the scanning strategy dependence of suppression of the power at the even multipoles in the case 
of the combination of the two surveys. Note that the power at odd multipoles is different from zero only for a ^ 0; this is of 
particular relevance in the case of £ = 1 (dipole), and may have implications also on data calibration with the dipolo and dipole 
modulation. See also the text. 



of the simulation) is absolutely negligible in practice in 
this context. 

Much more remarkable is the APS found when one 
of the two surveys is not complete. In practice, this may 
be the case in realistic experiments (even for the most 
symmetrical NSS) because of the different location of 
the various receivers on the focal plane. A particular 
care to this aspect should then be taken in the analysis 
of low multipoles. 

As analytically expected, for q 7^ a non-negligible 
dipole contribution also appears, as evident in the map 
in Figure El It has the same amplitude in the two sur- 
veys, as predicted. Note that in this case the APS is 
different from zero also for the odd multipoles, in agree- 
ment with the analytical results (see Figure Clearly, 



the effect is in any case very low, appearing only to 
higher order in a. 

For relevant values of a, a case that cannot be ex- 
haustively treated through the analytical approach, the 
power of the dipole and of the other odd multipoles sig- 
nificantly increases (see Figures El and . Note also in 
FigureElthe presence of remarkable features introduced 
by the significant angle between the directions of the 
main spillover and of the spacecraft spin axis. 

We note that an approximation for the scaling with 
a of the APS with £ — 1 can be derived by the analyt- 
ical expressions in equations 113.46^ and l|8.47|l . holding 
in the limit of small values a, replacing a with sin a 
which implies the sin^ a dependence evident in panels 
a2 and a3 of Figure|Hl Note also that this dipole scaling 
holds almost independently of the details of the scan- 
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ning strategy as clear from panels b2 and b3 of Figure |H| 
(see also the next subsection), because of its large an- 
gular scale origin. 

5.2.2 Cycloidal scanning strategies 

The complication of the observational strategy implied 
by a cycloidal option has a clear impact both on the 
final maps and on the APS. 

We first consider slow precessions (6 month period, 
Nc = 2). In the case a = 0°, the most relevant differ- 
ence with respect to the case of the NSS appears in the 
combination of the two surveys. As evident in Figure |H1 
structures appear on various angular scales because of 
the violation of the balance between the TOD of the 
two surveys implied by this observational strategy. The 
APS found in this case (see Figure (HJ is similar to that 
found in the case of the NSS but in the absence of com- 
pleteness of one of the two surveys. 

In the case a = 10° we find very similar results, 
except for the expected presence of the dipole and of 
the other odd multipole terms. 

Faster precessions (1 month period) with a signifi- 
cant value of Q (45°, see Figures |7| and (HJ , imply a re- 
markable dipole signature, analytically expected, a rel- 
atively more efficient smoothing of some large scale fea- 
tures in the combination of the two surveys, and then a 
relative decreasing of the power at multipoles less than 
~ 10 with respect to the previous case. 

Keeping a = 0°, we exploit the impact of different 
precession periods (see Figures |7] and |HJ . Note in Fig- 
ure|7|how the precession period appears in the maps as 
a "finger print" ; this effect is evident in general for Nc > 
2. This results is analogous to that found for the nu mber 
of hit s in scanning strategy analyses iPupac fc Taubeil 
l200fJ l ■ because different precession periods imply differ- 
ent periodic modulations in the dipole amplitude as ob- 
served by the main spillover. The comparison between 
panels bl, cl, c3 in Figure |H1 shows that, except for the 
quadrupole, the first weak and broad low £ bump in 
the APS (properly, in terms of Cgtil -I- l)/27r) migrates 
towards higher multipoles for increasing Nc- This is a 
direct effect of the decreasing of the angular scale of 
the "finger print" structures for decreasing precession 
periods. 

Finally, we have investigated on the impact of the 
precession angle amplitude by considering also an an- 
gle of 5° (panel c2 in Figure |HJ, instead 10°. We find a 
certain decreasing of the power at multipoles less than 
20-30 with respect to the previous cases, in particular 
jointly considering the two surveys. This can be under- 
stood on the basis of a simple continuity argument: de- 
creasing the amplitude, the cycloidal scanning strategy 
tends to the NSS. 



5.2.3 Symmetrical main spillover 

For sake of completeness, we have carried out also 
some simulations assuming a far sidelobe response az- 
imuthally symmetrical with respect to the direction of 
the main beam . We assume again a far sidelobe with 
a Gaussian profile, but in this case this refers only to 
the dependence on the colatitude from the beam centre 
axis. We assume again FWHMms = 20°, but it refers 
here only to azimuthal cuts. We renormalize the max- 
imum sidelobe response to have again a given ratio p 
(= 0.01 in the current computations) between the side- 
lobe response integrated over the whole solid angle and 
the dominant main beam integrated response. In the 
simulations, we consider two cases, with the maximum 
sidelobe response located at an angle of 90° from the 
main beam (i.e. a = 0°) or at an angle of 45° from the 
main beam (i.e. a = 45°). 

We considered again the NSS and cycloidal modu- 
lations of the spin axis, in this case with Nc — 12 and a 
semi-amplitude of 10°. 

As expected on the basis of an obvious symmetry 
argument, in the case a = 0° we find no effect indepen- 
dently of the adopted scanning strategy. 

For increasing a, the above full symmetry is bro- 
ken and in fact in the case a — 45° we find a non- 
vanishing effect. On the other hand, the azimuthal sym- 
metry adopted for the far sidelobes implies a strong sup- 
pression of the power at multipoles £ > 1 and the only 
surviving relevant term is the dipole as can easily 
understood by considering that for the case a — > 90° we 
recover the case of the original dipole pattern, smoothed 
and decreased of a factor p. The simulations carried 
out in these cases for a = 45° show a power at £ = 1 
very close to that found in the previous subsections for 
a = 45° for the same value of p. 



6 DISPLACEMENT OF THE MAIN 
SPILLOVER 

The relaxation of the assumption of main spillover cen- 
tre location on the plane defined by the spin axis and 
the telescope line of sight can be easily treated through 

This scheme is inspired by the case of observations taken 
with an antenna pattern similar to that of a feedhorn non- 
coupled to a telescope. 

In this case, the far sidelobe azimuthal symmetry, signif- 
icantly different from that considered in the previous cases, 
reflects into a suppression of the power at even multipoles 
and into a presence of a very weak power at odd multipoles, 
m any case smaller than ~ O.S/iK^ (and of ~ 0.01 - O-l^tK^ 
for ^ < 9) in terms of Ci£{£ + 1)/2tt, because of high order 

terms in a, analogously to the previous case (again the nu- 
merical estimates refer to p = 0.01). 



© 0000 RAS, MNRAS 000, 000-000 



16 C.Burigana, A.Gruppuso and F.Finelli 






Figure 9. Maps of DSC in the case of the NSS for small values of a and of the phase shift (top panels) and in the case 
of a cycloidal scanning strategy with small (middle panels) and large (bottom panels) values of a and of the phase shift. 
The parameters of each case are indicated above each map together with the reference to the panel in Figure Uni where the 
corresponding APS is displayed. Galactic coordinates and MoUweide projection are used. 



numerical simulations. Because of the spacecraft rota- 
tion about the spin axis, any displacement of the main 
spillover centre from that plane is geometrically equiv- 
alent to a proper phase shift (constant for all scan cir- 
cles) between the TOD containing the main beam centre 
pointing directions and the TOD containing the main 
spillover straylight signals. The TODs from the simu- 
lations described in the previous section can be then 
easily reordered to describe this case. 

Figure |5] shows the maps obtained (obviously as- 
suming a 7^ 0) for some representative choices of the 
relevant parameters while Figure [THl displays the corre- 
sponding APS (for instance, the phase shift refers here 
a "delay" of the main spillover straylight signals). 

Note that, as expected, in the case of small values of 
a and of the phase shift the results do not significantly 
change (see panels dl and d2 in Figure lTnt with respect 
to the case in which the main spillover centre is located 
on the plane defined by the spin axis and the telescope 
line of sight because the leading terms are still those at 
the lowest order in a and in the phase shift, as discussed 
above. On the contrary, for significant values of a and 
of the phase shift the odd multipole power (not only at 
= 1) is similar to that found for the contiguous even 



multipoles (see panel d3 in Figure \Wi for both a single 
survey and the average of the two surveys 

Analytically this phase shift can be parametrized 
generalizing equations 13.191 and 13.21)1 as follows 



9ms = COs{9mb T Z?) « + C , 

finis ='fis± sin(6'™6 =F /3) a -f O (q^ 



(6.1) 
(6.2) 



with P defining the displacement from 6mb and where 
upper (lower) signs refer to North (South) towards 
South (North) motion of the main beam. Of course, 
when /3 is vanishing, equations 13.191 and 13.201 are 
recovered. 

Replacing 16.11 and 16.21 in the definition of Isl 
[see equation 13.221 or 13.231 ] it is possible to compute 
the map 



Tf', 



(6.3) 



Here Tj,'^ is given by equation 13.331 [i.e. the introduc- 

In the case of small a and significant phase shift we nu- 
merically find again that the power at a given even multipole 
is larger than that at the contiguous odd multipoles (again 
except at £ = 1), although with a mitigation of its dominance, 
depending on the specific set of considered parameters. 
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Figure 10. APS of the DSC for various configurations. Panel 
dl refers to the NSS for a = 10° and a phase shift of 10° . 
The others panels refer to cycloidal scanning strategies with 
a semi-amplitude of 10° . Panel d2 refers to Nc = 2, a = 10° 
and a phase shift of 10°. Panel d3 refers to Nc = 12, a = 45° 
and a phase shift of 60°. In all panels triangles and crosses 
refer to a single survey and squares to the combination of the 
two surveys. Compare panel dl, d2, and d3 respectively with 
panel a2, b2, and b3 in Figure 13 See also Figure |51 Note the 
significant power at odd multipoles in the panel dS due to 
the relevant values of a and of the phase shift. Sec also the 
text. 



tion of P leaves uneffected the zeroth order, see equa- 
tions H6.1^ and Ij6.2|l ] while the first order is changed 
as 



tL (/3) = cos /3 tJ^ + sin (3 tL , 



(6.4) 



wiiere label " and ^ refer to tiie plane defined by tiie 
directions of tlie spin axis and main beam centre. The 
coefficient is given by equation H3.38|l and is 

tL = 2ci [-1 + i-ir] h - C23 [1 + (-1)"] h , (6.5) 
where C23 =02+03, and 



(6.6) 



dip{di cos ip — d2 sin ip)xYi^{6{x), ip) , (6.7) 



with 6[x) = arccos x. For symmetry reasons in equation 
16.511 only odd I contributions turn on. When m is odd 
tlien /2 = (because of tlie integration over x) and 7i 
can be non- vanishing (for I > 1, m has to be different 
from ±1); when m is even then 7i = (because of the 
integration over x) and I2 can be non-vanishing (for 
i > 1, m has to be different from 0). For ^ = 1, we 
compute 



and 



16 



2230(2 



(6.8) 



(6.9) 



We conclude that to first order in a only odd £ are 
present and the shift (3 modifies the dipole term and 
switch on ^ > 1 terms. When /3 ~ a than the shift is a 
second order term and to linear order does not appear. 
This is in agreement with the numerical result found in 
the linear regime in a, as evident from the comparison 
of panel a2 of Figure |H1 with panel dl of Figure \Wi (see 
also footnote 15). 

Finally, we give the average Tem over two or an even 
number of surveys in case of a shift P different from zero: 



Tem = a cos /3 Tl 



(6.10) 



Since, as already mentioned, 7^ only for £ = 1 (see 
equation 113.3811 and Subsection 13.81 see also Subsection 
13.71 , only the dipole term survives after two (or an even 
number of) surveys, also for /3 7^ 0. This result agrees 
with panel dl of Figure [Tni where /3 ~ q. 



7 IMPLICATIONS FOR SPINNING SPACE 
MISSIONS 

Numerical simulations confirm the analytical result that 
this effect is particularly remarkable in the case of an 
odd number of surveys, while a proper average of an 
even number of surveys greatly reduces its amplitude. 

On the other hand, for complex scanning strate- 
gies and or in the presence of uncompleteness of one 
of the considered surveys the suppression through aver- 
aging of this effect is significantly reduced. Clearly, the 
choice of a given scanning strategy is driven by many 
other aspects (mission constraints, sky coverage, redun- 
dancy, overall reduction of sy s tematic effects, and so on ; 
see e.g. iBernard et alJ j2002l) . lDupac fc TaubeJ J2005ft . 
iMaris e^dTl2o'o6h 'l. Our analysis shows that, in partic- 
ular at low multipoles, a special care should be taken in 
combining data from multiple surveys. Since, typically, 
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at low £ the sensitivity is not a problem in comparison 
to systematic effects and foreground removal, it could 
be preferable to avoid to include in the low £ analysis 
data exceeding an even number of complete surveys to 
reduce a priori the amplitude of this effect, or, at least, 
to compare the results found in this way with those de- 
rived by using the whole set of data. 

We have considered here the DSC in the case of a 
single receiver. This is the case of PLANCK-like optical 
configurations in which the sky signal is compared with 
a reference signal. In other optical configurations, for 
example in the case of WMAP, CMB anisotropies are 
mapped by comparing the sky signal in two different sky 
directions. Clearly, in this case the effective DSC comes 
from the difference between the dipole signal as seen by 
the far sidelobes corresponding to the considered pair of 
receivers. Therefore, the final effect will mainly depend 
on the different orientation and response of their main 
spillovers 

7.1 Calibration 

We note that for a 7^ 0, and according to the value 
of p, a non-negligible power at ^ = 1 may appear. It 
increases with a and can produce on the map peak 
values of ~ ±10/iK (for a of few tens of degrees) 
for p — 0.01. Clearly, the value is small compared to 
the kinematic dipole amplitude of ~ 3 mK, but, de- 
pending on p, it could be non-fuUy negligible when 
compared with the calibration accuracy of CMB space 
mission with both dipole and its ~ 0.3 mK modu- 



around the Sun 


Bersanelh et al.lll997HPiat et alJl2002l 


ICaDDellini et al.l 


20031: iJarosik et alJl2006ll. 



7.2 Removal during data analysis 

Clearly, provided that the far pattern is well known, 
it is possible to accurately simulate this effect for the 
mission effective scanning strategy or to apply deconvo- 
lution m ap making schemes w orking on the full beam 
pattern faarrison et alj|200d) in order to subtract the 
DSC from data, maps, and APS. The first strategy has 
been in fact adopted by the WMAP team to remove 
the overall (i.e. from the global signal sum of the var- 
ious components) straylight contamination for all the 
frequency channels during the data analysis of the three 
year data CJarosik ct al. 2006.) . 

Clearly, a precise subtraction of the DSC relies on 
a very accurate knowledge of the far beam. 

^® For example, if they simultaneously point to the same sky 
region with similar responses the final effect will be signif- 
icantly attenuated. On the contrary, if they simultaneously 
point to sky regions with opposite dipole signs, the effect will 
be amplified. 



For example, an error produced by an overall un- 
derestimation (or overestimation) of the main spillover 
response of 3 dB (or 2 dB, or 1 dB) is equivalent in this 
simple scheme to a ~ 100% (resp., 60% or 26%) error on 
the value of p (or, equivalently, oi fsL, Fsl, see equation 
14.241 1'). The residual contamination in the map and in 
\/CV can be obviously obtained by the results presented 
above rescaled by the same factor. Analogously, Fig- 
ures 13 and 2] can be used to understand the error on 
the quadrupole recovery. While for B < 0, even a poor 
knowledge of Fsl (i.e. of p) implies in any case at least a 
partial removal of this systematic effect (with respect to 
the case in which this subtraction is not applied to the 
data - obvioulsy, this subtraction has a meaning only 
for Ap/p < 100%), the treatment of the case _B > is 
more difficult because the final effect of this removal will 
depend on the true value of Fsl (or p), on the error in 
its knowledge, and on the intrinsic quadrupole Ci^^ , 
given the parabolic behaviour of y{FsL)- This calls for 
a particular care in the knowledge of the main spillover, 
at a level better than ~ 1 dB. 



7.3 Comparison with Galactic straylight 
contamination 



iBurigana et alJ lj20o3) presented an analysis of the 
Galactic straylight contamination at 100 GHz and 
30 GHz taking int o account the optical simulation by 
ISandri et al.l Jio^, in the context of Planck LFI op- 
timization work. 

For values of p ~ 0.001 (see the quantity /% 
in Table 3 in ISandri et ahl \2QQ4] ). at 100 GHz 
iBurigana et alJ (|2004l i (see § 5) found Ct£{£ + 
l)/27r ^ 10~^/iK^. Rescaling the values reported in Fig- 
ure |H] for p — 0.01 to the case of p ~ 0.001 we find 
Ci£{£ + 1)/2tv in the range ~ [3 - 0.03]/iK^ (or in the 
range ~ [0.1 — 0.01]/iK^) from lower to higher multipoles 
for the case of a single survey (or of the combination of 
two surveys). 

We then conclude that, at least at frequencies close 
to the minimum of Galactic foregrounds (70 GHz), the 
DSC is more relevant at low multipoles than the Galac- 
tic straylight contamination. 

Clearly, the different frequency scaling of CMB 
and Galactic foregrounds implies an opposite con- 
clusion at significantly lower and higher fre quencies 
ijDe Zotti et alJll999l: iBouchet fc GisDert /'l999^ . 

Finally, we observe that the APS of Galactic stray- 
light contamination in the case of Planck LFI also 
shows a power typically larger for the even mult ipoles 
than fo r the odd ones - see Figure 9 of Burigana et al.l 
ll2001al) -. analogously to what expected from our anal- 
ysis in the presence of a large dipole term in the con- 
sidered diffuse component for a far beam dominated by 
the main spillover feature not far from the spin axis. 
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7.4 Polarization measurements 

The CMB polarization anisotropy is typically observed 
by combining the signals from a set of polarization re- 
ceivers. For example, for differential radiometers as in 
Planck LFI or in WMAP the minimum set of receivers 
for polarization observations consists of four radiome- 
ters coupled to two feeds. Although the CMB dipole 
is not polarized, the DSC may affect polarization mea- 
surements because of differences in the intensity in the 
various receivers and, in the simple case of combination 
of data from four radiometers coupled to two feeds from 
the differences in the dipole straylight signals in each 
pair of radiometers associated to the sa me feed. There- 
fore, the method described in § 4 of iBurigana et alj 
l|2004 l can be applied here to provide simple numerical 
estimates. Clearly, also in this case the effective DSC 
will mainly depend on the different orientation and re- 
sponse of the main spillovers corresponding to each pair 
of radiometers associated to the same feed. Assuming a 
difference of a factor of ~ 2 between the corresponding 
signals in these two radiometers, we find that the con- 
tamination in the Q and U Stokes parameters is similar 
to that present in the intensity. 

This underlines the relevance of a very stringent 
control of this systematic effect for future accurate po- 
larization measurements, being the CMB polarization 
E mode at low multipoles less than few xlO~^/iK (in 
terms of {Ci£{£ + l)/27r)^''^) and being the expected B 
mode smaller at least of a factor of two or three, accord- 
ing to the c urrent WMAP uppe r limits on the tensor to 
scalar ratio dSoergel et alJl200d) . 



8 STATISTICAL ANALYSIS OF THE 
QUADRUPOLE 

The low multipoles of the CMB anisotropy pattern 
probes the largest scales of our universe, far beyond the 
present Hubble radius. The two experiments so far capa- 
ble to measure such low multipoles - COBE/DMR and 
WMAP - have detected an amplitude for the quadrupole 
surprisingly low compared to theoretical expectations. 
It is not clear if th e origin of this anomaly is sta- 
tistically signi ficant (lEfstathi ou 2003) , or related to 
foregrounds l|T ep^mark^ de Oliveira-Costa fc H amiltonI 
l2003l : lAbramo fc Sodrel2003l : ICopi et alJl2006D or possi- 
ble systematic effects like the one discussed here. In the 
following we shall discuss the impact of DSC on this 
issue. 

As discussed in Section |7| our modelling should be 
modified in order to investigate the DSC in the context 
of instruments based on differential measurements for 
CMB anisotropics, such as COBE/DMR and WMAP, 
which is beyond the scope of this work. However, we can 
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Figure 11. Distribution (in terms of counts; j/-axis) of the 
observed quadrupole C2 {a;-axis). 

make some considerations on the possible maximum 
amplitude of this effect at low multipoles as if measured 
by an instrument like the one we have modelled. 

First, we note that the value of the octupole has 
changed significantly from WMAP 1-yr release to the 
3-yr one, reducing the problem of unexpected low am- 
plitudes mainly to the quadrupole. Therefore, a system- 
atic as DSC which leaves unalterated the octupole is in 
principle in better shape after the WMAP 3-yr release. 

Second, as it can be seen from Section 4, the am- 
plitude and correlation of the DSC contamination can 
in principle reduce for the instrumental values chosen 
( which are realistic for Planck ) a theoretical expecta- 
tion of the quadrupole (C2 ~ 10"^ /iK^) to the one ob- 
served. Therefore, it is important to address the proba- 
bility of how DSC affects the observed quadrupole, since 
such a systematic can increase or decrease Cj'^'^ depend- 
ing on the sky realization (i.e. T.fjf^). 

We perform here a statistical analysis aimed at the 
computation of the pro bability of increase or decrease 
of the quadrupole (see ICruppuso. Burigana fc Finellil 
feOOdl for a preliminary analysis) . We have made 50000 
Gaussian distributed estractions of such that on 

average Cf ~ 10'^ /xK^ and we have computed the 
distribution for the observed quadrupole C2 (see Fig- 
ure ^^i- In Figure [T^ we show the distribution for the 
intrinsic quadrupole (i.e. Cf^^) and in Figure [T^ we 
plot the distribution for the pure DSC (i.e. the differ- 
ence C2 — Cf^^). In Figure IT^ we see that the dis- 
tribution of DSC has a Gaussian profile with a mean 
different from zero, as can be understood from equa- 
tions (14.221 and (14.231 1. We check that the found mean 
(202.3 fiK^) and standard deviation (404.3 ^K^) of this 
distribution are in agreement with the values expected 
respectively from equation 14.221 and equation (14.231 . 
for the chosen value of Fsl (= 4.59 x 10~^). The prob- 
ability of quadrupole amplitude increasing (decreasing) 
is 69% (31%). 

To this aim, we consider the case of an odd number of 
surveys. 
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Figure 12. Distribution (in terms of counts; j/-axis) of the 
intrinsic quadrupole C^^^ (a;-axis). 
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Figure 13. Distribution (in terms of counts; j/-axis) of the 
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pure DSC C2 - C|^^ (x-axis). 



9 CONCLUSION 

We have developed an analytical model and numer- 
ical analyses to evaluate the DSC in spinning CMB 
anisotropy missions. Although our study is mainly de- 
voted to the Planck project, the formalism and method 
are relatively general and allow to focus on the most rel- 
evant DSC implications. We quantify this systematic ef- 
fect as a function of few parameters: the relative power, 
p, entering the main spillover region with respect to 
the total one, the solid angle subtended by the main 
spillover region, the angle, a, between the directions of 
the main spillover and the spacecraft spin axis, and, 
finally, a phase, /3, describing the displacement of the 
main spillover centre direction from the plane defined 
by the spacecraft spin axis and the telescope line of 
sight. The first and third of these parameters turn out 
to be the most relevant (at least for small a, as in the 
case of Planck). Also, we have addressed the coupling 
of this effect with the observational strategy, with and 
without displacements of the spacecraft spin axis from 
the ecliptical plane. We have investigated the relevance 
of performing multiple surveys and the effect introduced 
by a possible uncompleteness (or overcompleteness) of 
one of the considered surveys. 

The analytical approach, applied to a simple ob- 
servational strategy (e.g. the NSS, in the context of 



Planck) and perturbatively in a, has been used to fo- 
cus on low multipoles, albeit some general properties 
have been found for odd and even i. The systematic 
effect vanishes for an even number of complete sky sur- 
veys, except for £ — 1 when a ^ 0. We predict a con- 
tamination of the dipole itself, independently of the con- 
sidered number of surveys and mainly depending on p 
and a. We have shown that when a phase /3 is intro- 
duced, not only the dipole but also all the odd mul- 
tipoles are present to linear order in a. We find also 
that the quadrupole is affected and its observed ampli- 
tude is related to the intrinsic sky realization {Ti^^). 
A statistical analysis has been performed aimed at the 
computation of the probability of increasing or decreas- 
ing of the observed quadrupole amplitude with respect 
the intrinsic one. 

With numerical simulations, based on a dedicated 
updated code already applied in straylight analyses, we 
verify the above results and extend the analysis to dif- 
ferent scanning strategies and significant displacements 
of the main spillover centre direction from the spin axis 
direction even for higher I and uncomplete (or overcom- 
plete) surveys. Various aspects relevant in CMB space 
projects (such as implications for calibration, impact 
on polarization measurements, accuracy requirement in 
the far beam knowledge for data analysis applications, 
scanning strategy dependence) have been discussed. 

Our analysis shows that DSC should be carefully 
taken into account for an high precision calibration of 
the data from the Planck receivers and for an accurate 
evaluation of the quadrupole amplitude of the CMB pat- 
tern. 
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